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Abstract 

We characterize the graphs formed by two-dimensional soap bubbles as being exactly the 3- 
regular bridgeless planar multigraphs. Our characterization combines a local characterization of 
soap bubble graphs in terms of the curvatures of arcs meeting at common vertices, a proof that 
this characterization remains invariant under Mobius transformations, an application of Mobius 
invariance to prove bridgelessness, and a Mobius-invariant power diagram of circles previously 
developed by the author for its applications in graph drawing. 
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Figure 1: Left: Boy Blowing Bubbles (1663), hyFians van Mieiis. Right: Girl Blowing Soap Bubbles 
(1674), by Pierre Mignard. Both images pubhc domain from Wikimedia Commons. 



1 Introduction 

Blowing soap bubbles has long been a popular pastime for children and a subject of art (Figure [T]), 
but the mathematics of soap bubbles has also been a subject of study and some mystery for many 
years, and the geometry and combinatorial structure of soap bubbles are not well understood. 
The Kelvin conjecture, stating that a foam with truncated-octahedron cells has the minimum 
total surface area among all soap bubble foams with equal-volume cells, was posed in the 19th 
century by Lord Kelvin and finally solved negatively in 1993 by the discovery of the Weaire- 
Phelan structure, a more parsimonious foam with two different cell shapes j35|[38]; proof that this 
newer structure is optimal remains elusive. A related problem, the minimum-area cap for the 



ends of hexagonal honeycomb cells, also remains open 15 . It is conjectured that minimum-area 
four-bubble clusters, for given amounts of gas in each bubble, form stereographic projections of 
a four-dimensional simplex, with each 2-face of the simplex projecting to a spherical surface in 
three-dimensional space, but this has not been proven, and there are clusters of six bubbles for 
which the numerically-computed surfaces separating the bubbles are not spherical [33] . Even as 
simple-sounding a problem as the double bubble conjecture (proving that the optimal structure for 
clusters of two bubbles has three spherical patches sharing a common circle) took many years to 
solve [17|[20] 



In this paper we study the special case of soap bubbles formed between glass plates, set close 
enough to each other that each bubble surface spans the entire gap between the plates. Effectively, 
these bubbles are confined to a two dimensional plane rather than three dimensional space. Such a 
system is familiar through the "soap bubble computer" , in which pegs are placed between the two 
plates at specified locations; a soap film that connects the pegs forms a minimal (though generally 
not minimum) Steiner network, providing a heuristic and approximate physical solution to an NP- 
complete problem [l||4|[Tl|[l9||2l]. However, we consider here bubbles supported only by the two 
glass plates, without additional pegs (in equilibrium, assuming perfectly impermeable bubbles, but 
not necessarily globally minimizing the total bubble length). When restricted in this way, bubble 
clusters may be described combinatorially by planar graphs, with vertices where three bubbles 
meet and edges for the interfaces between pairs of bubbles. With infinitely many bubbles of equal 
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areas, the optimal shape for the bubbles is a hexagonal tiling of the plane [15||16j, and similar 
tilings arise also for finitely many equal bubbles [g] , but other conditions give rise to more irregular 
arrangements of bubbles (Figure [2]) The question we consider is, with finitely many bubbles of 
variable areas, which planar graphs are possible? 

We provide a complete answer to this question, by 
showing that the graphs of planar soap bubble clusters 
are exactly the bridgeless three-regular planar graphs. To 
do this we use two sets of long-known facts about soap 
bubbles: Plateau's laws, describing the local geometry of 
the surfaces and junctions in a soap bubble cluster, and 
the Young-Laplace equation, which relates the curvature 
of soap bubble surfaces to the pressures of the gas within 
the bubbles. Based on these facts we prove two addi- 
tional properties of planar soap bubbles clusters: a char- 
acterization of them in terms of purely local conditions 
on the curvatures of the arcs meeting at each vertex, and 
a lemma stating that a Mobius transformation of a pla- 
nar soap bubble cluster is itself also a planar soap bubble 
cluster. The fact that a soap bubble cluster must form a 
bridgeless graph follows easily by combining this Mobius 
transformation property with the Young-Laplace equa- 
tion. We show that bridgelessness and 3-regularity are 
not just necessary but also sufficient by using two of the 
three steps in a recent method of the author for finding 




Figure 2: Soap bubbles between two 
glass plates. Public domain nega- 



tive image by Klaus-Dieter Keller from 
[Wikimedia commons! 



circular-arc drawings of planar graphs 12 . The first step applies to 3-vertex-connected graphs, 
and applies a novel Mobius-invariant power diagram of disks (defined using three-dimensional hy- 
perbolic geometry) to a family of disks constructed from the planar dual of the given graph by 
the Koebe-Andreev-Thurston circle packing theorem j32j. The second step generalizes from 3- 
vertex-connected graphs to 2-vertex-connected graphs by means of SPQR trees [7|. (The third step 
of this graph drawing algorithm, generalizing from 2-vertex-connected graphs to arbitrary graphs, 
does not work for soap bubbles.) Because of our use of power diagrams, our result can be seen 
as a validation of Sullivan's suggestion that "we might look for foams as relaxations of Voronoi 
decompositions" 
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Soap bubbles are closely related to Lombardi drawings of graphs, drawings in which all edges are 
represented as circular arcs meeting at equal angles at each vertex [8- 10, 12 , and as detailed above 
our characterization uses methods previously used for Lombardi drawing; however, not all planar 3- 
regular Lombardi drawings represent soap bubbles, as Figure [3] (left) shows. Our result can also be 
seen as an analogue for soap bubbles of Steinitz's theorem that the graphs of three-dimensional con- 
vex polyhedra are exactly the 3-vertex-connected planar graphs fsT] , and of several related theorems 
characterizing the combinatorial structure of geometric objects. Similar graph-theoretic character- 
izations (with different connectivity conditions) are also known for simply-connected polyhedra 



with axis-parallel edges 13 and for subdivisions of rectangles into smaller rectangles 23,24,36 



among other structures. Indeed, the Koebe-Andreev-Thurston theorem that we use here can also 
be viewed as a result in this vein. 
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Figure 3: Left: A collection of segments and arcs that obeys Plateau's laws, but does not form 
a planar soap bubble cluster, because it disobays the Young-Laplace equation (the curvatures do 
not add to zero at some of its vertices). Three of its segments do not separate different regions, 
impossible for soap bubbles. Modified from a figure in ^O]. Center: An irregular planar soap 
bubble cluster. Right: a soap bubble cluster with the topology of a buckyball. 



2 The classical conditions 

Two classical results on soap bubbles may be used to describe their geometry. Plateau's laws, 
observed experimentally in the 19th century by Joseph Plateau |28| and proven rigorously for 



minimal surfaces in 1976 by Jean Taylor 34 , state that in a three-dimensional soap bubble cluster, 

• Each two-dimensional surface has constant mean curvature. That is, the two principal cur- 
vatures of the surface take the same average value at all points in the surface. 



At each one-dimensional junction of surfaces, exactly three surfaces meet, and they form 
dihedral angles of 2tt/3 with each other. Such a junction is called a Plateau border. 



• At each endpoint of a Plateau border, exactly four Plateau borders and six two-dimensional 
surfaces meet, and the borders form angles of cos~^(— 1/3) with each other (the same angle 
that would be formed by two rays from the center of a regular tetrahedron through two of its 
vertices). 

The Young-Laplace equation, formulated in the 19th century by Thomas Young and Pierre-Simon 
Laplace, states that the mean curvature of a two-dimensional surface in a soap bubble cluster 
is proportional to the difference in pressure on the two sides of the surface, with a constant of 
proportionality determined by the surface tension of the fluid forming the soap bubbles [22j . 

In the case of a soap bubble cluster formed between two flat plates, the bubble walls are 
perpendicular to the plates, and all cross-sections of the bubbles by a plane parallel to one of 
the plates are congruent to each other. In this case, the principal curvature in the direction 
perpendicular to the plates is zero, and we may simplify both laws to describe the planar figure 
formed by the cross-sections: 

• The figure consists of one-dimensional curves of constant curvature; that is, circular arcs or 
line segments that do not cross each other. 

• At each endpoint of one of these arcs or segments, exactly three curves meet, and they form 
angles of 27r/3 with each other. 

• The curvature of any one of the circular arcs (the inverse of its radius) is proportional to the 
difference in pressure between the bubbles it separates. Bubbles with the same pressure as 
each other are separated by line segments, with zero curvature. 
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If a system of bubbles and bubble boundaries obeys these laws, the local forces on each vertex 
or arc caused by pressure and surface tension will sum to zero, so it will be in an equilibrium 
state (though possibly an unstable one [l4|[37]). Thus, we define a planar soap bubble cluster to 
be a finite collection of curves, obeying the planar restrictions of Plateau's laws and the Young- 
Laplace equation (for some assignment of pressures to each bubble). Equivalently, in graph-theoretic 
terminology, it is a planar embedding of a graph with circular-arc edges, forming angles of 27r/3 at 
each vertex, for which it is possible to find a pressure assignment to the bubbles that is consistent 
with all of the arc curvatures. As we describe in the next section, it is possible to replace this 
existential condition (the existence of a consistent pressure assignment) with a simpler statement 
about the curvatures of the three arcs meeting at each vertex. Examples of planar soap bubble 
cluster are shown in Figure [Sj center and right. 

3 Local characterization 

To state our local characterization of planar soap bubble clusters, we need a signed variation of 
curvature. Given a circular arc ending at a point p, of curvature c > 0, we define the signed 
curvature of the arc at p to be c whenever the arc curves clockwise as it leaves p, and to be — c 
whenever it curves counterclockwise. For line segments, of course, the signed curvature is zero. 

Lemma 1. A finite collection of circular arcs and line segments forms a planar soap bubble cluster 
if and only if it obeys Plateau 's laws and if, at each endpoint of an arc or segment, the sum of the 
signed curvatures of the three incoming curves is zero. 

Essentially the same result (for bubbles on a sphere instead of in the plane) can be found in 
Quinn Maurmann's appendix to [SO], so we defer the proof of the lemma to an appendix. 

4 Mobius invariance 

When three three-dimensional bubbles form a cluster in which all surface patches are spherical, it 
was known to Plateau that the three inner surface patches have centers of curvature that all lie on 
a single line |18|. Although centers of circles and collinearity of points are not Mobius-invariant 
properties, we can restate the same collinearity property for planar bubble clusters in a Mobius- 
invariant way, in terms of the crossing points of the circles. This restatement will allow us to prove 
that the Mobius transformation of a planar soap bubble cluster is itself another planar soap bubble 
cluster. 

Lemma 2. Let three circular arcs meet at angles of 211/2, at a point X, with centers of curvature 
Ci for i £ {1,2,3}, and let ri = \XCi\ be the radii of each of the three arcs. Then the following 
three conditions are equivalent: 

1. The sum of the three signed curvatures of the arcs is zero. 

2. The three center points Ci are collinear. 

3. The three circles with centers Ci and radii ri have two triple crossing points. 

Again, we defer the proof to an appendix. This triple crossing property is closely related to, 
and inspired by, the existence and uniqueness of the standard double bubble, for which see e.g. 
Proposition 14.1 of J27]. Figure |4] shows a set of of arcs, circles, centers, and radii meeting the 
conditions of the lemma. The lemma applies directly only to circular arcs and not to straight line 
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Figure 4: Arcs meeting the conditions of Lemma [2] 



segments, but the first and third conditions are equivalent more generaUy: For two arcs and one 
segment that meet at a point X, obeying Plateau's laws and having zero curvature sum, the two 
arcs must have equal and opposite curvatures and there are again two triple crossings between the 
two circles and the line containing the three given curves. It is not possible for one arc and two 
segments to have zero curvature sum, and the case of three line segments can be thought of as 
having a second triple crossing "at infinity" . 

The Mobius transformations are a family of transformations of the plane, augmented by a 
single point at infinity; if the points of the plane are represented by complex numbers, each such 
transformation is either a fractional linear transformation 

az + b 

z ^ 

cz + a 

(where a, 6, c, and d are complex numbers with ad — be ^ 0) or its complex conjugate. The 
Mobius transformation of a circle or line is necessarily another circle or line, and the Mobius 
transformations can be characterized as the largest group of transformations that preserves circles 
in this way. Mobius transformations are also conformal: if two curves meet at an angle 6, their 
transformed images also form the same angle. The Mobius transformation of a line segment or 
circular arc is generally another line segment or circular arc, but there are two additional undesired 
cases: it could instead be a ray, or a pair of oppositely-oriented rays on a common line. 

Mobius transformations do not preserve physically meaningful quantities of soap bubbles such 
as their pressures or the total length or area of soap film. They do not even preserve the finiteness 
of the bubbles: it is possible to find a transformation that takes the unbounded region of the plane 
to a bounded region and vice versa. Nevertheless, surprisingly, they transform soap bubble clusters 
into other soap bubble clusters: 

Lemma 3. Let B be a planar soap bubble cluster, and let t be a Mobius transformation with the 
property that no curve or arc of B is transformed into a ray or double ray. Then t{B) is also a 
planar soap bubble cluster. 

Proof. We first consider the case that neither B nor t{B) has any straight segments. By Lemma [l] 
and Lemma ^ B consists of arcs meeting at angles of 27r/3 such that each three arcs that meet 
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have circles forming two triple crossings. The meeting angles and the triple crossings of circles are 
preserved by Mobius transformation, so t{B) is also a collection of arcs meeting at angles of 27r/3 
such that each three arcs that meet have circles forming two triple crossings. By Lemma [2] again, 
the signed curvatures at each meeting point of t{B) add to zero, and by Lemma [l] again, t{B) is 
a planar soap bubble cluster. 

If B or t{B) contains straight segments, the equivalence between zero sums of signed curvature 
and triple crossing points (counting the point at infinity as a crossing of any two lines) at any 
endpoint of a segment can be shown even more easily, as detailed in the remarks following Lemma[2j 
and the result follows in the same way. □ 

Lemmas [2] and [3] are stated in different terminology as Lemmas 5.1 and 5.2 of ^39j. A hint of 
this Mobius transformation property of bubble clusters can also be found in Sullivan's conjecture 
that four-bubble clusters are stereographic projections of a 4-dimensional simplex [33], as any two 
such projections are related to each other by a three-dimensional Mobius transformation. 

5 Bridgelessness 

Our graph-theoretic characterization of planar soap bubble clusters asserts that the graphs of these 
clusters are bridgeless, or equivalently (since they are three-regular graphs) that they are two- 
vertex-connected. In the language of soap bubbles, this means that all of the arcs and segments in 
the cluster separate two different bubbles. As we now prove, this is necessarily true for a planar 
soap bubble cluster. 

Lemma 4. In a planar soap bubble cluster, it is not possible for a segment or arc to have the same 
bubble on both of its sides. 

Proof. It is not possible for a circular arc to have the same bubble on both sides, because the 
pressure would necessarily be the same on both sides and the circularity of the arc would violate 
the Young-Laplace equation. But then it is also not possible for a line segment to have the same 
bubble on both sides, because it is straightforward to find a Mobius transformation that transforms 
this segment into a curved arc, violating either the Young-Laplace equation or Lemma |3j □ 

6 A Mobius-invariant power diagram of disks 

In [12] we used three-dimensional hyperbolic geometry to devise, from a given set of disjoint disks 
(or complements of disks) in the plane, a partition of the plane into regions bounded by circular arcs 
that is invariant under Mobius transformations, and we used this partition to construct Lombardi 
drawings of 3-regular planar graphs. In the same paper we observed that the junction between 
regions for three mutually tangent disks could be found as the isodynamic point of the triangle 
formed by the three points of tangency, allowing the calculation of the diagram for this case to avoid 
hyperbolic geometry, but we did not find a purely two-dimensional description of this diagram for 
arbitrary disjoint disks and we did not describe how to extend it to disks that might intersect. Here 
we describe the same structure as a minimization diagram for a distance function defined between 
(point, disk) pairs in the plane, and we extend it to disks that may not necessarily be disjoint; that 
is, we show that it can be interpreted as a form of planar power diagram. Note that although this 
diagram is itself Mobius-invariant, the distance that it minimizes is not; this situation is much like 
that for the Delaunay triangulation of points on a sphere, which is Mobius-invariant even though 
spherical distance is not |3]. We keep our description brief and intuitive, for space reasons and 
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Figure 5: Left: The power distance from the red point outside the circle is the length of its tangent 
segment; the power distance from the blue point inside the circle is half the length of a chord 
bisected by the point. Right: The radial power distance from the red point outside the disk is the 
radius of the two congruent circles tangent to each other at the red point and both tangent to the 
circle; the radial power distance from the blue point inside the disk is the negation of the radius of 
the two congruent circles. 




because the additional generalization compared to [12] is not necessary for our results on soap 
bubbles. 

The three-dimensional hyperbolic construction from |12| 
is as follows. We consider the plane in which the disks lie 
to be the boundary plane of points "at infinity" for a three- 
dimensional Poincare halfspace model of hyperbolic space; the 
circle bounding each disk Di is the set of limit points of a plane 
Pi in this hyperbolic space. We form the three-dimensional 
Voronoi diagram whose sites are the planes Pi, and we extend 
this diagram from hyperbolic space to the boundary points 
of the model. The Voronoi diagram is invariant under con- 
gruences of hyperbolic space, and these congruences become 
Mobius transformations when restricted to the boundary of 
the model, so the resulting diagram is invariant under Mobius 
transformations as desired. 

To extend the hyperbolic construction to disks that may 
intersect, we consider each of the planes Pi in hyperbolic space 
to be the boundary of a halfspace Hi containing Di. We define 
the signed distance of a point q from Pi to be the positive 
distance from q to Pi if q is outside of Hi, and the negative 
distance if q is inside of Hi. Then, in hyperbolic space, the minimization diagram of these signed 
distances has bisectors that are hyperbolic planes, so each cell of the minimization diagram is 
an intersection of halfspaces. When the disks are disjoint, the previous construction (the Voronoi 
diagram of disjoint hyperbolic planes) is recovered as a special case. For a point q within hyperbolic 
space, its nearest neighbor by signed distance can be found (as in the Euclidean case) by considering 
concentric spheres centered at q and determining either the largest such sphere that is contained 
in a halfspace Hi or (if no such sphere exists) the smallest such sphere disjoint from all halfspaces. 
When g is a point at infinity of hyperbolic space, its nearest neighbor can be found in the same 
way, but using horo spheres in place of concentric spheres; these are shapes that in the halfspace 
model of hyperbolic space are modeled as spheres tangent at q to the plane at infinity. 

In the Euclidean plane, the classical power distance of a point g to a circle C is either the length 
of a tangent line segment from g to C (if q is outside C) or minus half of the length of a chord 



Figure 6: A packing of seven disks 
and one disk complement, with its 
Mobius-invariant power diagram. 



Modified from a figure in 12 
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Figure 7: Cross-section of a halfspace (blue) and 
a horosphere externally tangent to it (red) in the 
halfspace model of hyperbolic space. The plane 
on which the shadows fall is the plane at infin- 
ity of the model, and the black curves show the 
line qc (where q is the point of tangency of the 
horosphere to the plane at infinity and c is the 
center of the disk from which the halfspace was 
defined) and the two tangent circles in the plane 
at infinity that are congruent (as Euclidean fig- 
ures) to the cross-sectional circles. 

bisected by q (if q is inside C); see Figure [Sj left. If d is the Euclidean distance from q to the center 
of C and r is the radius of C, then by the Pythagorean theorem, in either case, the squared power 
distance multiplied by the sign of the power distance has the simple algebraic form d? — r'^. The 
power diagram is the minimization diagram of power distance (or of — r^); it has a polygonal 
cell for each of a given set of circles within which the power distance to that circle is less than the 
power distance to any other circle 

Analogously, in the Euclidean plane, given a point q and a disk D, define the radial power 
distance from q to D hy finding a pair of congruent circles Ci and C2, tangent to each other at q 
and both tangent to D, as shown in Figure [5j right. If q is outside D, the radial power distance 
is the radius of these circles; if q is inside D, it is the negation of the radius. Again, by the 
Pythagorean theorem, the radial power distance has the algebraic form (d^ — r'^)/2r. We claim 
that the Mobius-invariant diagram of disks defined above by three-dimensional hyperbolic geometry 
is the same thing as the minimization diagram for radial power distance. To see this, consider a 
planar cross-section of the three-dimensional hyperbolic space, defined by the plane in the half- 
space model that passes through q and the center c of D and is perpendicular to the boundary 
plane. This cross-section contains the closest point of D to q; the half-space bounded by D appears 
in the cross-section as a semicircle congruent to half of D, and the largest horosphere containing 
q and entirely inside or entirely outside of the half-space appears in the cross-section as a circle 
congruent to Ci and C2. Indeed, the cross-section of the half-space and the horosphere together 
may be formed from D and Ci or C2 by folding them upwards at right angles from the given 
Euclidean plane, along line qc (Figure [7|). From this it can be seen that the minimization diagram 
of radial power distance prioritizes the circles Ci in exactly the same way that the minimization 
diagram of hyperbolic signed distance prioritizes the corresponding horospheres, and therefore that 
the two diagrams coincide. 

7 Realizing graphs as soap bubbles 

We are now ready to describe the constructions from fT2] that we use here to generate a planar 
soap bubble cluster for any bridgeless 3-regular planar multigraph G. 

We consider three different cases, according to the connectivity and existence of multiple edges 
in the given graph G: 

• If G is 3-vertex-connected multigraph that is not a simple graph, then G must have two 
vertices and three parallel edges. For, any other connected 3-regular multigraph must have 
a pair of parallel edges between some pair of vertices u and v, each of u and v must also be 
adjacent to a singleton edge uu' and vv' , and the removal of u' and v' would separate u and 
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Figure 8: Left: double bubble representing a 3-edge multigraph. Center: gluing multiple Lombardi 
drawings together on an S node with alternating virtual and non-virtual edges (schematic view, not 
an actual planar soap bubble cluster), modified from a figure in 12 . Right: a planar soap bubble 
cluster that does not come from a circle packing. 



V from the rest of the graph, violating 3-connectivity. But a graph with three parallel edges 
is easily realized as a standard double bubble, with one straight line segment and two arcs 
spanning angles of 47r/3 (Figure [sj left). 

If G is 3-vertex-connected and has no multiple edges, its planar dual is a maximal planar 
graph G' . By the Koebe-Thurston-Andreev circle packing theorem [32] , we may find a set 
D of disks (or complements of disks) with disjoint interiors, corresponding one-for-one with 
the vertices in G', such that two vertices are adjacent in G' if and only if the corresponding 
two disks in D are tangent. By performing a Mobius transformation if necessary we may 
ensure without loss of generality that D includes one disk complement and that the other 
elements of D are disks. Let H be the Mobius-invariant power diagram of -D, as shown in 
Figure |6j None of the boundaries between regions in n can cross the disk complement in D 
and therefore no boundary can be a ray or a double ray. H has a region for each disk of G\ 



so by planar graph duality its vertices and edges form a drawing of G 12 . 

To show that H realizes G as a planar soap bubble cluster, let ABC be any triangle of 
the maximal planar dual graph G' , and let p, g, and r be the points of tangency of the 
corresponding three disks in D. An inversion ^ centered at the isodynamic point of triangle 
pqr (a special type of Mobius transformation) will take pqr to an equilateral triangle (this 
property characterizes the isodynamic point and is often used to define it) and therefore 
causes the transformed images of the three disks to be congruent. By the Mobius invariance 
of the power diagram, ^(H) is the Mobius-invariant power diagram of i{D), and in ^(H) (by 
symmetry) the three curves forming the boundaries between the regions for p, q, and r form 
straight line segments or rays that meet at angles of 27r/3. These boundary curves lie on 
three lines with have two triple crossings in the extended plane: one where the segments 
meet, and one at infinity. Back in H, the boundaries between the same three regions are 
the images under the transformation of these three line segments; therefore they are 
necessarily either line segments or arcs of circles, they again meet at angles of 27r/3, and they 
again lie on lines or circles that have two triple crossing points. By Lemma [2| they also have 
signed curvatures that sum to zero. We have shown that each three incident curves in 11 
satisfy the conditions of Lemma [T] and therefore by that lemma 11 is a planar soap bubble 
cluster realizing G. Figures [3] and [6] give examples of soap bubble clusters constructed for 
3-connected graphs in this way. 

Finally, suppose that G is bridgeless but not 3-connected. We may decompose G into 3- 
connected components, whose connections to each other are described by an SQPR tree 
(Figure loj) l7l[25]. In an SPQR tree, each node represents a 3-connected component of the 



9 




Figure 9: A planar 3-regular and bridgeless but not 3-connected graph (left) and its SPQR tree 
(right). The blue dashed curves link pairs of virtual edges. Figures from fTsj. 



graph. Each component is a smaller graph that may be a cycle (an S-node), a two- vertex 
multigraph (a P-node), or a 3-vertex-connected graph (an R-node). Within a component, 
some edges may be designated as "virtual" , and each edge of the SPQR tree connects two 
virtual edges in different components; the original graph G may be recovered by identifying 
and then deleting pairs of virtual edges. In the SPQR tree of a 3-regular planar graph, each 
P-node has three edges, each R node is itself 3-regular and planar, exactly one virtual edge 
of each linked pair must belong to an S-node, and each S-node must be an even length cycle 



that alternates between virtual and non- virtual edges 13 ,29 . 

To realize G as a planar soap bubble cluster, we separately realize each P-node and i?-node 
and then glue together the virtual edges of one S'-node at a time. To glue together the planar 
soap bubble clusters linked to an S'-node, we perform Mobius transformations of each cluster 
so that the circular arcs realizing their virtual edges lie on a common circle, and so that the 
remaining parts of the clusters are shrunken to lie within disjoint disks that meet the circle 
in the correct order. We then use the arcs from each cluster within the disk that contains it, 
and the arcs of the shared circle to connect them, as shown schematically in Figure [s] (center) . 
Details of this shrinking and gluing process may be found in [12|. Since each triple of incident 
curves in the resulting drawing is locally a Mobius transformation of a triple of incident edges 
in one of the 3-connected components, it again satisfies the conditions of Lemma [T] and again 
the whole drawing 11 is a planar soap bubble cluster. 

This gives us our main result: 

Theorem 5. A multigraph G can he realized as the arcs and junctions of a planar soap bubble 
cluster if and only if G is planar, 3-regular, and bridgeless. 

Proof. Planarity and 3-regularity follow from Plateau's laws. The requirement that a planar soap 
bubble cluster be bridgeless is Lemma |4j Conversely, the construction outlined in this section shows 
that every planar 3-regular bridgeless graph has a realization as a planar soap bubble cluster. □ 

This method can be made into a practical algorithm for realizing graphs as soap bubbles, which 



we have implemented in the 3-connected case 12 , but because it depends on the circle packing 



theorem it necessarily produces approximate numerical results rather than precise vertex and edge 



coordinates, and its running time depends on the precision of the results |5 26 . We note that, 
although our circle packing method is capable of realizing all graphs of planar soap bubble clusters, 
it is not capable of generating all valid geometries of planar soap bubble clusters. In particular, 
the clusters we generate have the property that within each bubble there is a circle tangent to all 
the arcs of the bubble, but this property does not hold for arbitrary soap bubble clusters; Figures 
[2] and [s] (right) provide counterexamples. 
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A Proofs of Lemmas 



Recall that Lemma [T] states that a finite collection of circular arcs and line segments forms a planar 
soap bubble cluster if and only if it obeys Plateau's laws and if, at each endpoint of an arc or 
segment, the sum of the signed curvatures of the three incoming curves is zero. Similar local 
characterizations were known already (see the appendix of [30^) but for self-containedness, we give 
here a proof. 

Proof. In one direction, in a planar soap bubble cluster, let x, y, and z be the pressures of three 
bubbles that meet at a point. Algebraically, {x — y) + {y — z) + {z — x) = 0, and therefore the three 
signed curvatures which by the Young-Laplace equation are proportional to x — y, y — z, and z — x 
also add to zero. 

In the other direction, suppose that finitely many circular arcs obey Plateau's laws and have 
signed curvatures adding to zero at each point where three arcs meet. We must assign pressures to 
the bubbles formed by the arcs, such that the curvatures obey the Young-Laplace equation. Assign 
zero pressure to the outside region. For each bounded bubble B, choose arbitrarily a curve c that 
starts in the outside region, ends within B, meets the arcs only at proper crossing points interior to 
an arc, and has only finitely many crossing points. The pressure within B can then be determined 
as the sum of the pressure differences determined by the curvatures of each crossed arc. 

To show that this system of pressures obeys the Young-Laplace equation, consider any two 
bubbles Bi and B2 separated by an arc A, and let ci and C2 be the curves from the outside 
region to Bi and B2 by which their pressures were determined. We may link ci and C2 by another 
curve that remains entirely within the outside region, forming a single (possibly self-crossing) curve 
that connects Bi to B2, and we may continuously deform this curve, crossing finitely many triple 
points of the system of arcs as we do, until it forms a line segment crossing arc A. Each time 
this deformation causes the curve to cross a triple point, the assumption that the curvatures at 
that point add to zero ensures that the sum of the pressure differences of crossings along the curve 
remain constant. Before the deformation, this sum was the pressure difference between Bi and B2 
in our system of pressures, and after the deformation, it is the pressure difference required between 
bubbles Bi and B2 for the curvature of A to obey the Young-Laplace equation. Since these two 
pressure differences are equal, A (and by the same argument every arc) obeys the equation. □ 

Lemma [2] states that the following three conditions on three circular arcs meeting at a point X, 
with centers of curvature Ci and radii (as shown in Figure [4]) are equivalent: 

1. The sum of the three signed curvatures of the arcs is zero. 

2. The three center points Ci are collinear. 

3. The three circles with centers Ci and radii rj have two triple crossing points. 

Similar statements can be found in the literature (e.g. see Lemma 5.1 of |39j) but again we supply 
a proof for completeness. 

Proof. We partition the proof into four implications between the three conditions of the lemma. 
(1) (2): 

Suppose that the three signed curvatures sum to zero; then one must have a different sign 
than the other two. Without loss of generality (by permuting the indices and by mirroring 
the configuration, if necessary) we may assume that the signed curvatures are 1/ri, — l/r2, 
and l/r^. By assumption, the sum of these three quantities is zero; by multiplying each term 
of this sum by rir2r3 and rearranging, we obtain the equation rir2 + r2rz = rir^. 
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Note that the three hues XCi form angles of 7r/3 with each other, bceause they are perpendic- 
ular to the arcs, which meet at angles of 27r/3. The assumption on the signs of the curvatures 
implies that angles C1XC2 and C2XC3 must both equal 7r/3, and angle C1XC3 must equal 
27r/3. Now consider the areas of the three triangles C1XC2, C2XC3, and C\XC^. In any 
triangle, the area can be computed by the side-angle-side formula as half the product of two 
adjacent side lengths with the sine of the angle formed by the same two sides. Thus, these tri- 
angle areas are |rir2sin|, ^r2r3sin|, and ^rirssin^ respectively. But sin7r/3 = sin27r/3, 
so the already-obtained equation rir^ + r-iVz = rir^ implies that triangles C1XC2 and C2XC3 
together have the same total area as triangle C1XC3. This could only happen if the three 
centers Ci, C2, and C3 are collinear, for otherwise the sum of the areas of C1XC2 and C2XC3 
would differ from the area of C1XC3 by the area of triangle C1C2C3, which is zero only when 
these three points are collinear. 

^(1): 

By (2) the three center points Q are collinear; assume without loss of generality that C2 
lies between Ci and C3 on their common line. Because the three centers Ci, C2, and C3 
form angles of 7r/3 rather than 2tt/3, the signed curvature of the middle center C2 has the 
opposite sign to the signed curvature of the other three curvatures; we may assume without 
loss of generality (by mirror reversing the configuration if necessary) that these three signed 
curvatures are 1/ri, — l/r2, and 1/^3. As in the previous case, the three lines XCi form angles 
of 7r/3 with each other so angles C1XC2 and C2XC3 must both equal tt/3, angle C1XC3 must 
equal 27r/3, and by collinearity the two triangles C1XC2 and C2XC3 together disjointly cover 
the same region of the plane as the single triangle C'lXC-s. 

We can apply the side- angle-side formula to this region in two different ways, giving the 
equation 

1 vr 1 TT 1 . 27r 

-rir2 sm - + -r2rs sm - = -nrs sm—. 

But the factors of 1/2 cancel, as do the factors of sin7r/3 = sin27r/3, leaving the simpler 
equation rir2 + r2rz = rir^,. Dividing all terms by rir2rz gives l/r3 -|- 1/ri = l/r2, and 
rearranging gives 1/ri — l/r2 -|- l/r3 = as desired. 

The three circles have at least one triple crossing point, at X. Let £ be the line through the 
three centers, assumed to exist by (2). Then because i passes through each circle center, a 
reflection across £ is a symmetry of each circle and therefore of the whole configuration of 
three circles. The point X cannot lie on ^, because two circles centered on a line cannot cross 
at a point that is also on the line, they can only meet at a point of tangency, violating the 
assumption that the arcs meet at angles of 27r/3. Therefore, the reflection of X across i is 
also a triple crossing point, and the three circles have two triple crossings as (3) states. 

^(2): 

For any two intersecting circles, it is necessarily true that their centers lie on the perpendicular 
bisector of the chord connecting their two intersection points. But by the assumption of (3) 

that there arc two triple crossing points, the chords defined by each pair of the three given 
circles coincide; therefore, their perpendicular bisectors also coincide in a line I that contains 
all three circles. 

□ 
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